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FOREWORD 


The vork dooumented In this report weu perfomed In the Advanced 
Space area of Vought Corporation for the RASA Johneon Space Center under 
Contract KAS9-1U776, Task 3.7 "Technology Assessment." The purpose of the 
vork has been to identify design requirements for future spacecraft heat re- 
jection systems and to evaluate the iiiq>aet of these requirements on the design 
of radiators. The study outlines development work needed to efficiently 
integrate heat rejection systems into future spacecraft having large heat re- 
jection capacities and long mission durations. Diq>hasis is given to conceptual 
designs which viU reduce the cost of the spEcecreft# 
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UTRODUCTIOi AID StMtARY 



Tb« objeetivat of bbii study art to Idsatlty doolpi roqulroMiits 
for futuro spsescrsft host rojoetloa systout to ovalusto tbo l«pset of tboso 
rsquirensnts on tba construotlon of eoovontlcnal puivsd fluid and bybiid boat 
plpe/pv^od fluid radiators t and to oonesptually dssi^i nsv boat rojeetion 
oystans vbieb nay inprove tbs psrfomanes cr rsdues tbs eoot of tbs opaeseraft. 
lbs study addrssBss bsat rsj action rsquirsnants vbieb ars largs eonparsd to 
those of existing systsns and nission durations vbieb ars relatively long. 

Heat rejection eapaeltles addressed in this study range fron 10 to 
250 KH. The operating life goal is arbitrarily selected as 3 years. The 
najor Inpaet of this on the construction of radlstor panels is that tubing vail 
thicknesses are sized so that there is only a snail probability of onteoroid 
penetration in a $ year nission. The radiator panel can be nade to survive 
longer nlsslons vith little vei|^ or cost inpaet by sinply increasing the 
vail thickness or by enqploylng neteorold bunpers. Other guidelines for de- 
sibling the heat rejection aysten are: transport fluid inlet tM^erature range 
fron lOO^F to 200^, transport fluid return tenperature range fron O^P to bO^P, 
and 200 Ib/Ktf electrical pover penalty for transport fluid puivlng. 

The study includes analyses of conventional punped fluid and hybrid 
heat pipe /punped fluid radiators vbieb detemine the optinun values of radiator 
fin thicknesses, spacing of heat pipes or transport tubes, Rsynolds nunbers, 
and other radiator design paraneters. Data neeessaxy for the selection of 
transport fluids and radiator surface coatings are also analysed in relation 
to future nission requirenents. In addition, nev radiator concepts are pre« 
sented vhlch are e:^eeted to be applicable to these requirenents. Oroundrules 
for deriving and evaluating the nev concepts give less eaq^tasis to traditional 
vei^t and deployed area than is curzvnt practice, and nore eq^iasls to other 
features such as stowage volune, user interface convenience, growth potential, 
adaptability to heat load varlationB, and reliability. 

The pazt of the study concerned vith the evolution of large heat 
rejection systems involves the specific conetzuetira of the elenents of the 
system; whether the radiators contain heat pipes, for exanple, as veil au 
procedures for constructing the system fron the elements. A large system 
built up from Independent subsystens with overslzlng, redundant conponents, 
and/or scheduled maintenance to compensate for element failures is signifl- 
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eantly lighter end nore feeeible taohnlcelly then one with e tingle treneport 
loop having the eaae cepecity and reliebillty« Dierefore* thie study evelu- 
etee elternetive ways to build up large qrateae tvcm aodulee by coaputing the 
eyeten else and probability of euoeaee taking account of the construction and 
reliability of the conponente. 

In euBBiary, the following approach le taken to iaprove future 
heat rejection eystens: 

1) elemente of the system are selected and desioned 
for optintum weight, cost, and reliability, 

2) new radiator coneepts are introduced that are 
uniquely desipied for future mission requirements, 

3) procedures for constructing large long-life systems 
are proposed that take advantage of modular con- 
struction, oversizing, and component redmdancy to 
ioq>rove performance and reliability* 


2.0 OFTINIZ/.TIOIf Of COMVIIFIGIIAL RADmOW 

Thi* ■•ota.on dltouini teilgn datallt for laprorlng oanT«iitlon«l 
■pao«erftft radlatori. It inolud«i data ralavant to ttoo ■•Itotlon of tranaport 
fluids and radiator coatings « rsllatlXltjr infomtlon on puip loop conponsnts* 
snd details of radiator panel construction for traditional parallel flow 
panels and for hybrid hsat pipe/pu«ped fluid pansls. 

2.1 PmiD SEUiCTlOM 


Technical Issues and property data for evaluating transport fluids 
in pusqped fluid or hybrid heat plpe/pmuped fluid radiator systems and working 
fluids in heat pipes are presented in this section. 

Transport fluids 

The selection of transport fluids is important since the fluid 


properties affect the radiator area, system weight, pressure drop, and com- 
ponent life. Of primary Interest are: (l) punqplng power penally (a function 
of viscosity, specific heat, at.d densl^), since this affects the puny> design 
and electrical power requirements, (2) heat transport eharacterlsties which 
determine the size of radiators and heat exchangers, (3) materials coa^ati- 
blllty and lubricity, which affect operating life, and (4) contamination 
threat. Also of inqportanee are availubility, freeze or pour point for low 
load or quiescent operation, flash point, toxicity, and flow stability. 

In previous studies pump power and conductance parameters have been 
developed ’ which are useful in assessing the relative merits of fluids for 
both laminar emd turbulent flow regimes. Ihese parameters were developed from 
basic fluid flow and heat transfer correlations by holding the boundary tem- 


perature, heat load, and system geometry constant. For these conditions, the 
puaqping parameter and conductance were foxmd to be dependent only on the 
transport fluid propez^les. Ihe pump power parameters fln*e: 


♦l 


P Cp 


for laminar flow 




p We fo** turbulent flow 
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and tita conductance pari 


itara ara: 


H for laninar flow 

L 

» §2 

n • JL a y for turtuXant flow 


vhera: |i • vlacoelty 

p ■ danaity 
Cp ■ apeciflc heat 
K ■ therntal conductivity 
p,. ■ Prandtl numher 

OMpufr of th. porfomac p«r»«te« of typlo,a cool»t. or. co»».d 

in Figurea 1 and 2. 

Fluida with low conductance parameters and high pump power para- 
■Otero ore poorly oultod for opococrdft cooling oyote» oppUcotlono. end 
vould bo ooloctod only If other flnldo felled boceuoo of other coneldoretione. 
The field eelected ueuelly ie not etmlght-foreerd beceuee flulde »lth de- 
.lrrt.le pawing power perenetere (euoh .. Freon 21) «ey heve lower oonductenoe 
pu^eutere then flnldo with poor paw pereMtere (ouch no glycol weter). Aleo, 
flulde which operete efficiently In the turbulent regime eaet be coweiud with 
aulde which fewor lemlner flow. Therefore, to eweluete the edepteblllty of 
cendldete flulde to e pertlculer epececruft It le ueueUy neceeeety to glee 
detelled coneldeiutlon to the redlator construction end the weight Increeee 
reiulred with the eelectlon of flulde with poor thermsl conductence, ee well 
es the liwect on the electrl.el power eyetem and the raw SoofO* 
selection of fluids with high pump power parameters. 

Properties of fluids which are candidates for future mission ap- 
plications are listed In Trt)le I. to optlssaa radiator design wUl exist for 
each nuld Which nlnlmlzes the system weight Including the weight of the power 
source for driving the pus®. To Illustrate how the choice of the tronsiort 
fluid affects a typical system design, calculations were mad. to deteimlne 
optimum system configurations for several candidate fluids assuming a weight 
penalty factor of 2l*5 Kg/KH. typical values of pump motor efficiencies 
end an assumed radiator environment of 0«F. Th. results, presented In Table 
II, Show that several fluids c<u> be eonsldered without seriously Impacting the 
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TABLE I 

PROPERTIES OF CAHDTIIato yrnrfifi 
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TABLE 11 

En-’ECT OF TRAM8P0RT FLUID ON WEIGHT AND GURFACE 
AREA FDR A KW SYSTEM 


TRANSPORT FLUtn 

MU20N 01 
FREON 11 
freon E-1 

FmON E-0 

KC-8S 

FC-7i> 

Fc-rr 

COOLANOL 1J> 

ETHYLENE GYLWL/WATER 
(RG-89a) 

OHONITE FC-lOO 


WEIGHT DELTA 
RELATIVE TO 
R-Ol (KG) 

AREA RATIO 
RELATIVE TO 
R-a (X) 

0 

1.000 

1.3 

1.00b 

0.? 

1.018 

2.1 

l.li»0 

0.5 

1.018 

0.0 

1.0b8 

3.1 

1.155 

1.0 

1.080 

- 0.0 

1.010 
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weight or ge«netrle configuration of the vehicle, and that the fluid selection 
can be based on other considerations such as naterials coaopatibility, operat- 
ing tenqperature range, vapor pressure, lubricity, and vehicle cootaaiinatlon 
in case of lecdiage. 

A property that is extremely isqportant in applications to pumped 
fluid or hybrid heat pipe/pui^ped fluid radiators is the fluid's susceptibility 
to flow instabilities. This property is a function of how the fluid viscosity 
varies with tesperature. 

In improperly desiffied spacecraft cooling systems flow Insttbili- 
ties may originate in the radiator panel where large changes in the transport 
fluid ten^rature ocexu: at low heat load conditions. For serpentine flow 
radiators with no parallel flow passages the flow insttbUity causes the pump 
to stall, and is usually accompanied by subsequent freesing of the transport 
fluid. This may damage the pump or motor, and in some orbit configurations, 
the astern may be permanently disabled because the transport fluid will not 
re-thaw irtien the operating condlticms change. 

Space radiators with parallel fluid transport passages asre also 
susceptible to flow instabilities at low heat load conditiems. Wh«a the dif- 
ference between the fluid inlet and outlet temperatures is too large, the flow 
distribution in the parallel tubes changes abruptly from being uniform, where 
each tube carries approximately the same flow, to non-uniform, where the flow 
in one or more of the parallel radiator tubes coBq>letely stagnates. Plow in- 
stabilities are undesirable when acccaapanied by freezing of the fluid in the 
non-flowing tubes since this may prohibit the flow from returning to the 
uniform distribution when the operating conditions change. Also, metal tubing 
may eventually fall due to stresses caused by the freezing and thawing of the 
transport fluid. To prevent flow instabilities the heat rejection system 
should be equipped with heaters or other methods of control which prevent un- 
favorable operating conditions from occurring, and a transport fluid should 
be selected which is resistent to flow instabilities. 

The properties of the transport fluid determine the range of 
operating conditions for which stable flow occurs. An equation from Ref. (2) 
gives the approximate operating limits for stable flow in terms of the trans- 
port fluid properties. The criteria for stable flow is: 
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v(T^) < 2 


i. 


vCp I 1 ♦ 

J 


irKIiu 


] dT 


r . t: 



( 31 ) 


Wiere: 

T • fluid outlet teiqperature 
o 

■ fl\ild inlet teiqperature 
T« ■ environment teiqperature 

V ■ fluid viscosity 

k ■ fluid thermal conductivity 
Cp ■ fluid specific heat 

V ■ radiator fin vidth 

e ■ radiator fin emittance 
o • Stefan Boltzmann constant 
Ru " Russelt Number 

Figure 3 gives the minimum outlet teiqperature for stable flow computed from 
Equation (l) for some cemdidate TCS transport fluids. The figure shows that 
RS-89a (ethylene glycol-water) has the more restricted operating range. This 
has been a significant factor which has limited the usage of ethylene glycol- 
water in space radiators. It does not eliminate the fluid from consideration 
for all applications because the thermal control system can often be designed 
to operate in the stable regime. However* the selection of a fluid such as 
Freon 21 would allow more freedom in the design of the control system. 

Heat Pipe Fluids 

The first consideration in the selection of a suitable wui’king 
fluid for heat pipes is vapor temperature range. A variety of characteristics 
must be examined in order to determine the most acceptable fluid for a proposed 
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FIGURE 3 


APPROXIMATE STABILITY CURVES FOR CANDIDATE RADIATOR 


FLUIDS 





•PPlicutlon* Sow of tho prim ro^ulromnts arot 

1) Plttid/natarlai conpatlbillty 

2) fluid thermal at ability 

3) Wetability of wick -and wall materlale 
Thermal properties (vapor preseuret 
latent heat* conductivity, viscosity , 
and surface tension) 

A convenient mans for quickly comparing wosking fluids is pro- 
vided by the liquid transport factor defined as: 

. • '■i. 

"I 

■ liquid transport factor 

* liquid surface tension 
“ latent heat of vaporization 

• liquid density 

■ liquid viscosity 

Curves of versus temperature are available in the literature for a nunflwr 
of candidate fluids. Data for candidate fluids for the teaq>erature range 
studied here are given in Figure k. At any selected operating temperature the 
fluid which exhibits the highest liquid transport factor will generally yield 
the best heat transfer characteristics. 

Contamination Threat of Fluids 

Evaluation of the contamination threat bf a heat pipe or transport 
fluid involves several considerations. One is the condensation temperature 
and energy associated with desorption or evaporation. This would be a factor 
in considerations of deposition rates on cold surfaces. Another could be the 
fluid's infrared absorption spectrum, which would influence the tolerable 
contamination of infrared sensors in som applications. A third ia the poten- 
tial of the fluid vapor for chemical reaction with other spacecraft materials. 

Table III gives materials compatibility data for several candidate fluids, and 
is taken from Reference 19 . 

Selection of Fluids for Radiato r Design Optimization stniHwa 

In following sections, detailed radiator optimization analyses are 
performed which require fluid property data. Freon 21 is the coolant employed 


N, 


where: 


N 
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FIGUKE 4 


HEAT PIPE FLUID TRANSPORT PARAMETER 








TABLE III MATERIALS COMPATIBILITY S 
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for the transport fluid In these studies. It's thermal properties result in 
minimum system weight, and Its operating twi^ierature range Is broader U»an most 
other fluids. It has good chemical stability, and low nonvolatile content. For 
heat pipes, ammonia is baselined as the working fluid. It has the highest liquid 
transport factor of the fluids applicable to the temperature range of this study. 
Other trade study criteria such as chemical stability and material compatibility 
appear to be acceptable, 
ii . THEHMAL tX)NTROii COATIMOB 

Thermal control coatings are an Important issue in this study because 
of the required operating life of 5-10 years, and the emphasis is placed on iilgh 
radiator performance and minimum cost. Good thermal control coatings sliould have 
the following properties: low solar absorptance (for sun viewing surfaces), high 
emittance, low outgassing, stability in space environment, and avoidance of static 
charge buildup (high altitude orbits). The discussion presented in this section 
Is a summary of approinrlate prior studies, and is based mainly on Reference 19. 

Candidates considered include: silver and aluminum backed Teflon, 
optical solai' reflectors, astroquartv. fabric, front surface coated Kat>t.on film, 
aluminum oxide coated aluminum, zinc orthotitonate paint, and a NASA-Goddard in- 
organic conductive yellow paint. Tnltial values of the emlssivity and sola* 
absorptivity, and degraded values of solar absorptivity are given in '.I'nble IV 
for some of the coatings. 


Ullver backed Teflon has become a very popular thermal conhrt)l t;oat- 
Ing because of low absorptivity (abou'. 0.1), high emlssivity (about O.O) and 
pm-ported stability In the space environment. Numerous spacecraft and i)rnbes 
have flown with this material. Figure 5 presents test and flight data ft>r silver 
Teflon showliig an initial rapid degradation (probably due to contamltrvUon ) fol* 
Jewed by minimal further increase in solar absorptivity in low earth tublf. For 
high altiiude fliglits where significant fluxes of electrons and preti.ic; are pre- 
sent in combination with solar ultraviolet, degradation continues at n t«te Llmt 
would appear miacceptable for a five year mission even If an asymptt-'l l •: tailing 
off eot?iirs. Heseurch for the Air F’orce has shown that silver Teflon lotien H.h 
mechanical Integrity and degrades such that Its solar absorptivity Increases by 
a factor of about i for a slmul.ated 5 year exiosure to high altitude ••i bits, There 
is also an Indication that low earth orbll degradation (solar ultraviolet only) 
mai Mint- be iiune severe for a 5 year miss h.n that would be ant I c I pnt« d . w. ii! 
with asLro>jnur'./, luaed silica fabric c<.iatlngs ^ Indicates tlml v;lien in 


























conjunction with a sllvcir Teflon aubatrate a raaeonably low aolar abaorptlvlty 
of about 0«lB la obtained, and Improved environmental atablllty la expected. 
Although the uae of aatroquartz la not an eatabllahed technology and It haa the 
dlaadvantage of being difficult to clean, allver Teflon with aatroq^uartz la a 
candidate for both low and high orblta becauae of Ita unique promlae of long 
life. Bare allver Teflon ia probably aui table for low orblta. Additional in- 
formation and performance data for metallized Teflon and for aatroquartz fabric 
are given in Tablea V and VI reapectlvely. 

Another candidate coating for direct aolar expoaure at high altitude 
la the Optical Solar Reflector (OSR) described In Table VIZ, consisting of silver 
or aluminum backed quartz tiles about one Inch square. Because these tiles are 
expensive, heavy, and must be bonded to the radiator they eu^e basically unde- 
sirable and not suited for large radiator systems. However, data indicates they 
have the lowest solar absorptivity available and are stable in the combined ul- 
traviolet, proton and electron environments of deep space. Thus these coatings 
should be considered in concepts which entail small radiators with significant 
solar exposvure. Their contamination threat from the adhesive must be evaluated. 

No paint coatings (organic or Inorganic) are know which cure stable 
In high altitude orbits. A new coating, zinc orthotltanate, is under develop- 
ment but insufficient data are available to assess its applicability at 
this time. An inorganic yellow paint Is also being developed by Qoddeurd for 
hi£^ altitudes. The only other known potenticdly stable hi^ altitude coatings 
aire evaporated silicon or aluminum oxides over metallized substrates. Table 
VIII gives data for these coatings. Figure 6 compares solar absorptivity data 
for candidate high altitude coatings. 

An alteimate approach idilch offers greatly Improved confidence is 
to orient the radiator such that very little solar irreuiiatlon is present. 

This would not only retard degradatlcm but would minimize. If not eliminate, 
degradation effects. In this case a siaqple durable coating such as anodized 
aluminum could be used (to provide a high emlssivlty). An added advantage 
would be the elimination of any contamination threat from either the coating 
itself or cm adhesive securing it. 

Data obtained during the past few years indicates that film sur- 
faces, such as multlloyer insulation or silver Teflon, can obtain a large 
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100 A OVERCOAT OF INDIUM/TIN OXIDE INDICATED TO RELIEVE 
MINIMAL SOLiiR ABSORPTANCE EFFECTS; Aa ■= C.02 (INITIAL) 
DEGRADATION OF OVERCOAT Aa = O.Oi IN 3000 ESH 
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MINIMAL SOOHCE OP OUTGASSING PRODUCTS 
DISADVANTAGES 

DIFFICULT TO AVOID OR REMOVE GROUND CONTAMINATION 
- APPARENI- GREATEB SETOITIVm TO OfiBITAl CMWAMIBATlai 
than METAtUZEO SILICA TILES OR TEPLM) 
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static charge buildup during solar substorm activity with high altitude space- 
craft. A static discharge can occur which damages the coating and provides 
spurious signals to the spacecraft. A significant amoxint of research has been 
directed toward alleviating static buildup in coating materials. A potential 
inherent advantage of the flexible radiator fin materiea also being developed 
under this contract is its potential to dissipate static charge. Since -ttie fin 
consists of a fine silver wire mesh embedded in the Teflon film, it is c<mceivable 
that this fin material, perhaps covered with astroquartz for stability, would offer 
a lightweight advanced technology approach to obtain a satisfactory high altitude 
radiator. Other countermeasures which may be taken to avoid static charge build- 
up include: 

e Addition of transparent conductive film to coating 
surface (induim/tin oxide, others) 

e Metal grid applied to svirface 

e Inh^ently conductive coating (astroquartz, conduc- 
tive paint) 

Selection of Thermal Control Coatings 

Ihe following reccsnmendations are made for selectii:g radiator coat- 
ings for long diuration missions: 

1) Low Eai*th Orbits — Metallized Teflon is the best currently 
qualified choice based on Initial a, e, and degraded properties. Contamination 
of the coating has the most adverse effect on o degradation. The radiators 
should be oriented away from thrusters and other sotarces of contamination when 
possible. 

2) ttigh Earth Orbits - No satisfactory coating available, althovigh 
astroqiaartz is probably acceptable with design allowances for degradation. 
Laboratory data indicate metallized Teflon is unacceptable without an overcoat 
to avoid static charging, although some satellite flight data suggest the pro- 
blem may not be as severe as would be expected. Optical solar reflectors have 

- -les4iu-swer-e-degradatlon than Teflon, but wei^t and cost make them unattractive 
for large areas. AstroqtJiartz fabric and Qoddard conductive yellow Inorganic 
paint have high initial values of solar absorptivity, are difficult to clean, 
and would be expected to degrade to a {W 0.3 in 5 years. More data is needed 
for other coatings such as the flexible radiator fin material, front surface 
Kapton film, zinc orthotitanate paint and alvimlnum oxide coated alianinum. 


2h 


2.3 


COMPOHEIIT LIFE STUDY 


The objectives of the component life study were to establish the 
current state-of-the-art conqponent lifetimes, to identify component designs 
which are capable of 5~10 yeco* operating lives, and to identify components whicli 
will require fiirther development. The approach taken wm to sur> ■ , '.he long- 
life cooq?onent technology* Much of the information presented in this section 
was obtained from Reference 19* 

Table IX summarizes operating life data obtained for typical space- 
craft coolant system components. The table shows that, with the exception of 
thermed control coatings, all components have projected life times of 5 years 
or greater. However, except for the fluid accmulator, five year lifetimes 
have not been demonstrated* Therefore, based on projected component perfor- 
mance and design edlowances for coating degradation, 5 yeeur operating life is 
achievable. The probability of success for 5-10 year missions depends on how 
the system is designed; whether it has been oversized or contains redundant 
elements for example. The details of how to design thermal control systems is 
discussed in Section 5*0. This section addresses only the operating character- 
istics of components of the system. 

Pumps 

Conteu:t8 with the major manufactui'ers of aerospace pumps indicate 
that the development of pximp motor assemblies capable of operating continuously 
for more than five years will challenge the state-of-the-art. Pump motor 
assemblies applicable to spacecraft thermal control systems have been tested 
for more than 20,000 hours without failure, and several designs have projected 
operating lifetimes greater than 5 years (Table X). 

Figure 7 paresents a typical pump selection curve which indicates 
regions of best efficiency for gear, vane and centrifugal pumps. The curve 
shows that centrifugal pun^s favor low pressure rise and hi^ flowrates. Often 
environmental control systems reqxdre a higher pressure rise, which could be 
met more efficiently by a gear or vane pump. However, gear and vane pumps 
generate contaminants and have more moving parts in contact and are thus un- 
desirable from life considerations. In addition they create pressure pulsations 
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and vltoatlons* To avoid these problems and still obtain satisfactory effi- 
ciencies (on the order of 30 %) double ended centrifugal pumps have been devised* 
These halve the pressure rise per stage enhancing overall efficiency* The 
double ended feature not only increases efficiency bub also unloads axial 
thrust Increasing lifetime potential* Further gains in lifetime may be possi- 
ble through the use of carbide bearings* vlth a good prospect for achieving 
the 3 year (44*000 hours) desired lifetime* 

Pump Motors 

Submerged canned stator AC motors have been proven reliable on space 
programs such as Apollo* Biosatellite* Shuttle Orbit er* and Saturn over a number 
of yeeurs* The technique of enclosing the motor stator in a sealed stainless 
steel can has been fully developed* to the point that the difference of effi- 
ciency between this and other designs is usually unmeasurable « The submerged 
motor has the inherent advantage of direct heat transfer to the fluid for 
efficient dissipation of motor heat* continual lubrication of the bearings* 
and the ability for hermatic sealing of the unit if necessary* Since the 
possibility of motor failure is generally associated with stalled punqps or 
motor overheat* and since the possibility of stalling a centrifugal ptinp is 
very remote* and the likelihood of motor overheat is also remote with direct 
fluid cooling* the integral submerged motor concept appears to be a leading 
loxig life choice* Solid state DC to AC inverters facilitate the use of 
canned stator AC motors on satellites. 

Accum\ilators 

The functions of the accumulator are to provide system pressure 
adequate to prev^t transport flviid vaporization at any teoqperature encountered, 
to provide positive punq> suctioa pressure* and to acconanodate fluid volume changes 
due to temperature excursions and leedcage. Two types of accumulators might be 
considered: bladderless accumulators and metal bellows. Two pressurization 
techniques* a gas pre-charge system and a regulated accumulator temperature 
system might be considered* The gas pre-charge system will require a larger 
overall volxime but will not require power. Vendor Information on accumulators 
is given in Table XI. 

Coolant Loop Bypass Valves 

The life of the coolant loop bypass valves is not considered to 
be a significant problem area* Because of the relatively steady heat load 
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typical of .pacocri>ft. tho «lv« »1U probatljr haw to cyola leaa than 50.000 
tlnaa In a 5 year Maalon. Spae. 4 nallfl«l alaotronlcaaly controlled yalwa haw 
haan toatad by Von*t for tha Bpaca Bhnttle program without difficulty for com- 
parahla duratlona over 30,000 cyolea. The daalgn Ufa of the Shuttle valw on 
tha Flow Control Aa.ambly (Fa) 1. 20.000 cyclaa. Blamanta of tha wlw which 
night fall include natal haUowa, teflon aaata, and etappar notora. Tha hallowa 
haw been axtanalwly taatad pact 50,000 cyclaa wharcaa Teflon aaata haw bean 
taatad to 500.000 cyclaa. Sypaaa valw atappar motora are capable of more th« 
30,000,000 atapa In a vacuun envlrcnnant. Thernoatatlcally actuated valvaa w 
wry low failure rataa are alao under davalopnent for auch prograna aa tha 

Satellite Infrared Experiment (SIRE) a 
Fluid Svivels 

The operating life depends strongly on the external loading an e 
aalactlon of aaala M ccag«naata for fluid tan®arat«ra varlatlona. With proper 
daalgna to account for aide and aUal loading, an operating Ufa of flw ya«a 
ahould be poaalbla. However, It la difficult to provide ayatan redundancy for 
rotating fluid awlwla bacauaa It la nacaaaary for all of the awlwla to be 

located on the axis of rotation a 
Heat Pine Ufe 

The majority of pUbllahed literature on heat pipe Ufa tecta Indl- 
^taa that Uf.tlmaa on tha order of 5 y«m^= «« 

of a heat pipe that la capable of contlnuoua operation over a 5-year period 

reuulrea metl.uloua attentlcm to aewral critical datalla. Some of 

heat pipe deaign detalla which can affect uaeful Ufetime Include: fluid chemical 

atabUity, fluld/materlal caapatlblllty, fluid purity, and metal cleanltaeaa. 

control during heat pipe manufacture alao cannot be owr-eng,haalaed in 

its relationship to continuous heat pipe performance. 

J,, the abow Itema m.»t be addreaead In order to Ineur. agalnat 

the primary cauaea of faUur. - corroaion and generation of non-condanaabl. 
gaaea. If the waU or wlch material la aoluble In the working 
tranafer la likely to occur between the condenaar ».d ewporator. Solid 
material wlU be depoalted In the ewporator reaulting In local hot apo a or 
blocking of the porea In the wick. Hon-condenaebl. gaa generation la probably 
the moat common cauae of heat pip. failure aa the non-cond«..ablea accumulat 

in the condenaer aectlon blocking fluid tranafer. 

Other lawortant life conalderationa In high perfor«H.ce heat pipe, 
include aurvival of the wick when aUbJected to the launch vlbroacouatlc apectrum. 
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Fatigue ilfe of the wick la Iwportant in dealgne requiring flexing* auch aa 
a girtbal aeetion. The topic of accelerated life teating and extrapolation to 
predict long tern performance la alao algi. ifloant to long life heat pipe 
aeaeaenanta. 


2.4 MICRQMBTBOROID PROTBCTKM 


When dealgnlng themal control ayatena for long duration nia alone* 
nicrometeorold conalderatlona have a algalfioant impact cc the dealgn of the 
radiatora and the coolant loop linea exposed to the apace environment. In 
heat pipe radiator dealffia the transport fluid manifolds must he shielded, 
and the radiator panel must he oversiaed to account for the loss of heat pipes 
hy meteoroid puncture. In pumped fluid radiators the manifolds and parallel 
flow cross tubes must he designed to withstand mlcrometeorold Intacta. Typi* 
cally the exposed cross section of the transport loop Is about the same for 
punned fluid or hybrid heat pipe/pumped fluid designs* so that the mass of 
armor or meteoroid bumper required is essentially independent of the radiator 
design. 

The transport loop wall thickness must be sufficient to retain the 
transport fluid pressure after being struck by a mlcrometeorold. The depth of 
the crater left by the most damaging meteoroid expected to strlice> the tubing 
during the designed operating life of the radiator is ccmputed from a ballis- 
tic equation which is based on ground test data* and a meteoroid flux model 
derived from penetrations of metal foils in near earth orbits. 

The mechanics of hypervelocity impacts are not conqpletely under- 
stood. An eqviation given in Reference (6) correlates the existing data 
reasonably well. The equation is 


t 


, 1/8 p 1/2 

0.65 if-) (rS- ) (V J 

®t h ® 


7/8(^ )19/10 


( 2 ) 


where: 


m 


thickness of material penetrated (cm) 
percentage elongation of teurget material 
mass density of sheet materleO. (gn/cc) 
mass density of meteoroid (gm/cc) 
nomal Impact velocity (km/sec) 
meteoroid diameter (cm) 
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Equation 2 generally predicts greater depths of penetration than other eqtia- 
(7 9) 

tions^'*'' derived for penetration of metals « end viU therefore provide for 
a conservative thermal control system design when ai^lied In conjunction with 
standard meteoroid environment models. 

Ihe radiator panel tubing vail thickness needed for stirvlval of 
a given s^vtem design d^ends on the projected target area, the exposure time, 
and the required probability of success for the mission. The probability of 
no meteoroid penetrations is given by 


P 


o 


^-CAtN 


(3) 


where: C 

A 
T 
N 


shielding factor (% of area e^^sed) 
projected area (m^) 
time of exposure (sec) 
meteoroid flux for particles capable 
of penetrating tubing (Particles/m^««ec) 


The tube vail thickness is computed from Equation (2) so that the 
number of particles determined from the meteoroid flux model having sufficient 
energy to penetrate the tubing gives the desired probability of success fraa 
Equation 3. 

For a five year mission the tubing must be designed so that only 

relatively large meteoroids are capable of penetration. For large meteoroids 

(9) 

the cumulative meteoroid flux model for sporadic and stream meteoroids is ^ 


logj^gN = - 1U.37 - 1.213 log^^m ik) 

Vhere N is the flux density for meteoroids having mass greater than or equal 
to ”m". The meteoroid mass may be expressed in texms of quantities in Equa« 
tion 2 as follows: 

“ ■ ? 

Equations 2 • ^ may be solved to give the tubing vail thicknesses required for 
a given design and required probability of success. 

Because of the danger of possible contaminaticm by leaking thennal 
control system fluids, the tube vail thickness should be selected so that there 
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is only a very small probability of a single meteoroid pvmeture. 

To minimize the contamination threat emphasis should be placed on 
designing the system to minimize the possibility of a single meteoroid penetra- 
tion with the least system weight and coiq>lexity. Table XII compares meteoroid 
penetration characteristics data for two typical 1?0 ft^ radiator designs as 
calculated from Eq,uations 2-5* The results show that high probabilities of 
mission success are possible with moderate radiator tube wall thicknesses* 

A significant weight savings is possible if plastic meteoroid 
barriers are employed so that the wall thickness of the transport tubing or 
heat pipes may be reduced. The elongation tern in Equation 2 is much larger 
for plastics (e « 300) than for metals (e « 3)* Limited experimental data in 
the literature substantiates the importance of this tern. 

Qualitative data for plexiglas and polycaibonate^^®^ shows that 
polycarbonate (p » 1*2, e « 5) for retenrding meteoroids. Hyperveloci'^ test 
results for polyethylene (p « 0.9, e » 500) in Bef. (7) are predicted conser- 
vatively by Equation (2). The equation predicts slightly greater depths of 
penetration than are measured at velocities below 12 Km/ sec, and much greater 
depths of penetration than Ref. (7) indicates will occur at velocities typi- 
cal of micrometeoroids, viz. 20 Km/ sec. 

Tests ”ere conducted at Texas A&M University^ in support of this 
contract to substantiate that Equation 2 gives the correct depth of penetration 
for FEP Teflon tubing. The results given in Flgiire 8 show that the penetration 
depths of nylon projectiles, p » 1.15, are predicted adequately by Equation 2 
whereas steel projectiles, p «* 8.0, penetrate greater depths than the eq'uatlon 
predicts. This indicates that the exponent of the projectile density in Equa- 
tion 2 may be too small. If this is the case the equation is conservative for 
predicting damage by micrometeoroids (p « 0.5)* 

Table XII shows that the weight of the radiator panel may be re- 
duced by more than 20 Ibm if Teflon meteoroid barriers are used. An alterna- 
tive procedure for reducing weight involves the application of meteoroid 
bumpers. Multiple wedl, or bumper type meteoroid protection systems are known 
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Radiator Tube Spacing 
Tubing I.S. 

Radiator Surface Area (Each Side) 
(2-Sided) 

Tubing Target Area 
Probability of no Penetration 
Tubing Wall Thickness Required 
Tubing Wei^t 

Wei{^t vith Teflon Meteoroid Barrier 


6 ” 

0.1875” 

150 Pt2 

1»».73 Pt2 
0.99 
0.139” 
51.3 Ibm 
29.6 Ibm 


12 ” 
0.50” 
150 Pt 2 

19.64 Pt2 
0.99 
0.151” 

55.6 Ibm 
31.4 Ibm 
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to have elgnlfloant might advantage over single wall or armor type ayetemst 
hut eonplioate the radiator design such that they have not been engdoyed in 
short missiMi duration applications where the mass associated with meteoroid 
protection is not significant* 

Ground test data have shown that the optimum multiple wall systan 
contains a thin first wall called a bumper, and a second or main wall >dilch 
is spaced at a sufficient distance from the bumper to permit a meteoroid 
striking the bumper to be dispersed before making contact with the main wall* 
Systems with more than two sheets ^ not provide as much protection as those 
with two properly designed sheets^ Honeycomb or other filler materials 
betweax the walls is less effective than a vacuum or air gap* The thickness 
of the bumper should be su:^iclent to cause damage to the meteoroid, but should 
not be so thick that material broken loose from the bumper has enough mass 
to damage the main wall. The optimum thickness of the bximper is approximately 
one-tenth the diameter of the projectile. 

Several esperical equations have been proposed to represait the 
relationship of the required thicknesses of the first and second walls, the 
spacing between the walls, the target material properties, and the projectile 
properties. A general form of the equation for the thickness of the main wall 
is 


K V ®S^p 

SI SI 


^d« 

m 


( 6 ) 


vhere K, o, g, y» and 6 are enpirlced constants, and S is the spacing between 
the plates. Table XIII lists the values of the exponents in Eq. 6 from various 
correlations in the literature. The table shows that there cure significant 


TiffilE XIII EXPONENTS IN EQUATION 6 


BEFEBENCE NO* 

o 

g 

Y 

a 

(tl/dm) 

11 

1 

- 2 

2 

k 

? 

12 

0 

- 2 

1 

3 

0.1 

13 

1 

- .5 

0.5 

1.05 

0.3 to 0.5 

Ik 

.270 

- 1.39 

7 

2.92 

0.25 to 0.1» 

15(a) 

? 

- 0*75 

0*6 

1.0 

0.25 

15(b) 

? 

- 1.1* 

0.6 

1.0 

0.10 
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dlff«reno 0 t In ibt form of tbo propoiod equntloni* Part of tho dlicrtpuoy 
is dut to vnrltttlons in tho thioknessos of the huaper sheoto vhloh hove o 
large Influenoe on the tjrpe of diautge to the main vail* Hhe eq,uations from 
Ref. (12) and Ref. (l^-h) are for optimum front sheet thiokhesses» and there- 
fore should not he aj^lied to more general oases. Of the remaining equations » 
the one from Ref. (ih) gives the most oonservative estimate of the material 
thiokness required to prevent penetration of the main vail* Ihls equation 
also predicts depths of penetration masured in Refs* (15) and (IT)* there- 
fore, Sqs* 7 and 8 helov, from Ref. (Ih), are ei^ployed in this vork to design 
radiators vlth meteoroid humpwrs. 


1/2 „O.270. 2.92 
t« - 5.33 (-^ ) ^ \ 


(7) 


O.ld^ < \ < O.k dju 


( 8 ) 


the materials constant, K, vhich vas originally measured for 
202l^-T3 aluminum targets and pyrex-glass projectiles, has heen modified to 
account for density variations of the target and projectile* 

When applying Bq* 7 to design radiators, there are tvo variables, 
tg and 8, vhich may be adjusted to provide meteoroid protection, idwreas in 
single vail relator designs, there is only one variable, the tubing vail 
thickness. The thickness of the bungper layer is fixed by Eq* 6 and the dia- 
meter of the smallest meteoroid that vlU penetrate the double vail system. 
Since there are tvo variables defined in Equation 6, one of the variables can 
be arbitrarily selected to minimise radiator fabrication difficulties, or to 
reduce vei^t. The radiator optimization computer routine described in 
Section 2.5 is programmed to minimise vei^t by solving Eq. 8 simultaneously 
vlth 

• 0 (9) 

non-linear expressions are obtained for tg and S idiloh must be solved itera- 
tively. 
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2*5 OPTIMIZATIOW OF PUMPED FLUID SPACE RADIATORS 

Manifold CoaflguratlonB 

Typically • punned fluid space radiators are designed with parallel 
flow passages which are Joined together by manifolds at the edges of the radia- 
tor. This design has been shown to be effective for ECS radiators with heat 
rejection capacities of the order of 10 KW. For large radiators, the system 
welj^t can be reduced If the flow routing Is modified so that the lengths of 
the manifolds and crosstube are not constrained by the overall radiator 
dimensions. 

For the configuration shown In Figure 9, the physical d^ tn emlo ng 
of the cross tubes are not tied directly to the radiator size so that the 
designers have more freedom to select the tube diameter, the Beynolds number, 
and other radiator parameters when optimizing heat transfer, pressure drop, 
and es^tosed ctoss section to the mlcrometeox^id environment. In effect, 
an arbiti*arily large radiator is constructed from a number of optimum sized, 
parallel connected panels having common manifolds. 

Very large radiator systems might re<iulre a more general panel 
configuration with parallel headers. In this case the manifold are 

independait of the radiator size, and may be varied to obtain an optimum design. 
This generality is not re<iulz«d for the radiators considered In this study. 
Statistical analyses presented in Section U.O show that very large radiator 
systems should be constructed from smaller subsystems having Independent 
pumped fluid loops. 

The parallel manifold configuration of Figure 9 Is Justifiable for 
the radiator sizes considered herein. The optimum cross tube diameters are 
relatively small because of micrometeoroid considerations so that close mani- 
fold spacing is necessary to minimize fluid pressure losses. Figure 10 gives 
the percentage radiator weight reduction that occurs with multiple manifolds 
for typical ECS operating conditions. 

The minimum radiator weight in Figure 10 occurs when the peuiel h gg 
four manifolds. Itois is typice^ of all the systems analyzed in this work. 

The wei^t of the headers becomes large in radiators with more than four mani- 
folds, and effects the advantage of the smaller tube diameters which are 
possible with additional manifolds. The radiator weight is relatively insen- 
*l'tlve to the number of manifolds for heat rejection rates of the order of 
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FIGURE 9 

PARALLEL MANIFOLD RADIATOR CONFIGURATION 



FREON 21 TRANSPORT FLUID 
P(0) ■ 0.99, 5 YEAR MISSION 
Tj^ - 100-F, T^yy « 40«F, T 




FIGURE 10 

WEIGHT REDUCTION FOR TYPICAL PUMPED FLUID 
SPACECRAFT RADIATORS WITH MULTIPLE MANIFOLDS 


10 Kw» so tbat the panel design ean be biassed vltbovt Height penalty by 
oonstralnts suoh as eould ooeur when adapting the panel to an existing 
deployawnt neohanlsa for exasqile* 

lolsation Caleulations 


Table XIV lists the design variables for the radiator panel of 
FlguM 9* To detennine the best radiator designr one nnst select an optimi- 
zation criterion in the fora of sone function of thei.dasign variables vhich 
is to be nininlzed, establish the relationship betveen the dependent and in- 
dependent variebleSf and then detemlne the values of the independent variebles 
^ch alnlaize the criterion function* in this section a procedure is des- 
cribed for ninlnizlng a function of radiator panel Height and area 


P • Ht + A (10) 

The Height of the radiator is considered to be the sum of the Heights of the 
manifolds, headers, cross tubes, radiator fin and coating, radiator fluid, and 
vel^t penalty for pressure drop* 

*t . W„ ♦ Wg ♦ ♦ Wp (U) 

other associated radiator system coiQtonenb Heights ean be included in the 
analysis by modifying constants in Eqs* 10 and 11* For example, the Height of 
the deployment system is proportional to the radiator area, and can be accounted 
for by increasing the constant Cx in Eq. 10. 

The manifolds and headers must be sized to give a uniform flow 
distribution. This is done by requiring that the pressure drop in the headers 
and manifolds be small in comparison to the pressure drop in the cross tubes* 

- OjiPjj (.12) 

IPg ■ CjlPg (13) 

For turbulent flow in the manifolds and cross tubes Eq. 12 becomes 
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The sum In Eqt* l4 nay be approxlaeted ae follows 
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Equation l4 beoones 
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The exterior aahlfolds carry only half as aiueh flow as the Interior unlfolds, 
and therefore have snaller dlaaeters 


(EXTERIOR) 
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0.772 d 


(17) 


For laminar flow In the manifolds and cross tubes, 

«2 0*25 

®ct ®2se 

For turbulent flow In the manifold and laminar flow In the cross tubes 
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The vel^t of the manifolds per unit area Is 
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( 20 ) 


The diameter of the hec^r la the same as the diameter of the manifold for 
designs with up to four manifolds. When there are more-than four manifolds. 
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The veigbt of the headers per unit area is 



" ITO ^ - “t/A) 

The vel^t of the cross tubes per unit area la 
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The vpJd thickness required to prevent failure by microDeteoroid puncture Is 
coaq>uted from the meteoroid environment model, a ballistic equation for depth 
of penetration of the tubing material, and the required probability of sur- 
vival. For long duration missions the tubing wall is sufficienUy thick that 

only relatively large meteoroids tore capable of penetration. The environment 

(9 ) 

model for large meteoroids is 


- 1.213 log^^M 


(2l») 


Where N is the average number of micrometeoroids with mass greater than on 
granes which will strike 1 m? of radiator area per second. Equation 2h assumes 
an average meteoroid density of 0.^ gm/cc. Thus the meteoroid diameter is 


djj » 1.7578 • 10"^ 1,-0 *2740 


(25) 


,( 6 ) 


A ballistic equation for depth of penetration of the radiator tubing is 




19/18 


t a 6.32 

The probability of no meteoroid penetration is 
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( 26 ) 


P(0) a e" (27) 

Where (A is the total exposed cross section and t is the time of exposure. 

The relative thicknesses of the cross tubes, manifolds, and headert. are deter- 
mined by minimizing tte radiator weight for constant survivability. The weight 
of the flow passages is 


® ^CT^CT * Vm ■*' Vh 


( 20 ) 


From Sqia. (29)* (26), and (27)» ttaa proloAlilllty of no penatratlon la 



For minimum weight 

d(«t) . ♦ A^dt„ ♦ Y^H • ® <30) 

For eonatent aurvlvahUlty, 
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Multiplying Ba* 31 hy an arbitrary constant (a Lagzmnglaa multiplier) and sub- 
tracting Sq. 30 
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For Eq,. 32 to be satisfied for all values of the multiplier 
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!Ehus the minimum weight occurs when the cross tubes, manifolds, and headers 
have equal thicknesses. The required thickness is 



(34) 


Where the exposed cross section Is related to the radiator area by 
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(35) 


similar equations derived for the main vail thlchness and bumper spacing in 
meteoroid bunq)er designs are programmed Into the radiator optimization computer 
routine listed in Appendix A. The radiator area Is computed by summing the 
radiating areas of each cross tube 
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(37) 


and R is the resistance to heat transfer betveen the transport fluid and the 
base of the radiating fln» If the resistance to conduction through the tube 
vail Is negligible* 


1 

vKfNu 


(38) 


Combining Eqs. 36 and 37 
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substituting Q = ^2*^^S ^®in ~ '^out^* integrating 
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The veight of the radiator fin and coating per unit area is 
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a. tMctae.. Of th. rodlrtor fin 1. » f»,ctlon of the «dl.tor fin nfflcloncf. 

For fin .fflcl«ol.. ttnu. O.T. tk. fining oquitlon roproNnko tht 

correlation obtained numerically ty Lieblein'^®', 


t ■ — 3 0 c 

KCl-n)r.8307 ♦ .T925(l-ti) ♦ .T050(l-n*))* 

The coetlng thlckn... 1. . nation of the t„e of eontlng (e.g. eU«r bneked 
Teflon) <md the enleeletty of the enrf«». m thl. rtndy, ellwr becked Teflon 
coetlng mth en e«i«.ltlt, of 0.8 1. eeecMd. tte coetlng tUckne«i, 
adheeivet for a two sided radiator is 0e012 inche 
The fluid weight per unit area is 
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Kheie 4P 1. conputed from friction fector wren. Beynold. number dete, «>d 0, 
is the pun 5 > power penalty factor (lb/KW)e 

Equation 44 is the final equation needed to determine the value of 
the function of weight and area that is to be minimised to obtain the optimum 
radiator desigUe To determine the best radiator design, trial values of the 
Independent variables in Table XIV are substituted into the equations listed 
above until the minimum value of the function P is established. Since this 
Involves many calculations and a logical search through the ranges of several 
Independent variables, a computerized approach is necessary. Listings of 
the radiator optimization routines developed for meteoroid armor and bumper 
protected desi^ie are listed in Appendices A and B. 

OFTIMIZATIOW of HIBRID heat PIPE/FUMPED PmiD RADTatobr 

Manifold Configuratlcma 

The manifold design of hybrid heat plpe/pua«)ed fluid radiators 
must take into account watt-lnch limitations of heat pipes, teii 5 ,erature drop 


49 


acroas the Interface from the transport fluid to the heat pipe vorking fluid, 
pressure drop In the transport flxild loop, and mlcrcmneteoroid protection re- 
quirements. Transport fluid flow routings such as are shown in Figure 11 
provide large radiating stirface consistent with watt-inch limits of heat pipes 
and minimum exposed manifold cross sections to micrometeoroids. The overall 
width and length of the radiator can be varied by changing the diameters of 
the heat pipes to adjust the edlowable distance between manifolds, and by 
selecting an Integral number of manifolds consistent with the approximate 
overall dimensions required. 

To minimize the temperature drop from the transport fluid to the 
heat pipe, heat conducting fins such as shown in Figure 12 may be employed. 

In this design ccm^act heat exchanger core is enq>loyed to transfer heat from 
the transport fluid to conductive fins attached to the heat pipes. The fins 
have a dual fvinction as heat conductors and meteoroid armor for the transport 
fluid loop. 

The manifold detail in Figure 12 is one of several analyzed in 
this work. It is competitive in performance for scane heat rejection require- 
ments with other designs studied, and would be relatively easy to fabricate. 
Other designs developed \mder Vought fimding are considered proprietary, 
and, therefore, are not discussed in detail here. Performance characteristics 
of some of the proprietary designs are ccanpared with those of the design of 
Figure 12 and of pumped fluid radiators in Section U.O. 

Heat Pipe Designs 

Exotic heat pipe designs are not required or desirable for hybrid 
reuliator applications. Since the cost of the heat pipes is a significant part 
of the total cost of the radiator, and high watt-inch capabilities are not 
required, the simplest possible heat pipe designs should be selected. Axial 
grooved ammonia heat pipes have acceptable thermal performance characteristics 
for most applications, cure highly reliable, and ccui be fabricated at less cost 
than other types of heat pipes. Center core wick heat pipes have slightly 
better thermal performance and cost slightly more than grooved pipes. Artery 
type heat pipes sliould be avoided because of their poor reliability and hlgli 
cost. 

Since each heat pipe in a hybrid heat pipe/pumped fluid radiator 
operates at a different tenq^erature , there is a slight performance advantage 
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FIGURE 12 TYPICAL HEAT PIPE RADIATOR MANIFOLD DETAIL 
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FIGURE 13 VARIABLE LENGTH HEAT PIPE RADIATOR 



possible if the lengths of the heat pipes are adjusted in accordance with their 
heat rejection requireiaents. Ihe problem to be solved in optimizing che design 
is to determine the length of the heat pipes in Figure 13 as a function in the 
flow direction (or fluid temperatures). The solution obtained by calculus of 
variations in Appendix C is 

t = ct’3''^ 

yh_e^e C is a constant. Equation shows that for nl^ 1 v^^nm n radiator area, the 
lengths of heat pipes near the transport fluid entrance to the radiator should 
be shorter than those near the outlet of the radiator. However, c€U.cuiations 
in Appendix C show that the area of a radiator with optimum heat pipe lengths 
is less than O.ljf smaller than the 8u:«a of a radiator with the same heat 
rejection capacity having constant length heat pipes. Thus, since vai-iable 
length heat pipes complicate the radiator design, and do not significantly 
improve performance, only constant length heat pipes are considered in this 
work. Computer analyses performed under a separate contract show that variable 
spacing of heat pipes also complicates the radiator design with less than 0.1? 

increase in performance. Thus only constant spacing geometries are considered 
in this work. 

The diameters of heat pipes must be sized so that they ore capable 
of transporting heat at the maximum reqvilred rate. For min-timim wei^t the 
heat pipe diameters should be as small as possible. Therefore, since the heat 
pipes near the transport fluid entrance are hotter, and reject more heat tlian 
those near the outlet, their optimvua diameters should be larger. However, 
since variable heat pipe diameters also complicate the problem of fabricating 
the radiator panel, only constant diameter designs are considered here. Jf 
future designs should become so weight critical that variable heat pipe dia- 
meters would be considered, it is probable that an alternative approach oi 
varying length or spacing in accordance with watt-inch requii-ements would be 
more attractive from a manufacturing viewpoint. 

Hybrid Heat Pipe Radiator Optimization 

The optimisation of hybrid heat pipe radiators is olmilur to thel 
of pumped fluid radiators in that values of independent variables such tu 
heat pipe spacing, and radiator fin thicknesses ore determined to minimize h 


i 
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funetion of voi^t and area* An additional oon^ieatlon la involved in nodel- 
Ing the heat transfer in the nanifold* Fbr exaiq>let the themal awdel of the 
nanil'old in Figure 1? must account for heat transfer in the eo^pact heat ex- 
changer core* the meteoroid shleldt the Joint between the manifold and the 
heat pipe* and the evaporative surface of the heat pipe* Pressure drop in 
tbs co8q;>aet heat exchanger core must also he accounted for* A general themal 
model is set up for each type of manifold design which includes adjustable 
parameters such as the heat exchanger core width "s" in Figure 12* In the 
optimisation analyses the best values of the adjustable parameters are deter- 
mined by computing the weight and area of the radiator for a range of permls- 
sable values of the parameters* The best manifold parameters and radiator fin 
parameters are detemined simultaneously to obtain a true optimum design* 
Oomptater routines were developed under Vought funding to 
optimise hybrid heat pipe radiator designs* These routines are considered to 
be proprietary* and are not discussed in detail here* Results of the optimi- 
sation emcdyses are presented In Section U.O vhere candidate radiatoxw based 
on puuqped fluid and hybrid heat pipe designs are eonqpared. 
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3.0 


ADVANCED RADIATOR C0NCEFT8 


In this section advsnced jradlator concepts are presented that are 
expected to he iq^pUcahle to the requlrenents of future spacecraft. These 
concepts give enphasls to adaptability, reliability, end growth potential for 
reduced development costsi and consider traditional wel^t and radiating area 
to be of secondary Importance. Oroundrules, desirable features, eusd techni- 
cal Issues for evaluating new concepts are listed In Tables XV, XVI, and 
XVII. Exanple applications of the advanced heat rejection systems are: the 
imall power module, the laxsa power modide, and experiments for space proces- 
sing, communications, solar power satellite demonstration, Ilfs sciences, and 
space construction demonstration. 

3.1 Modular. Self Contained. Long Life Radiator Astern 

The central Idea of the modular, self-contained, long life radiator 
system shown In Figure l4 Is to provide a modular deployable system with 
easily added or replaced modules. The heat rejection capacity can be easily 
ejgpanded as requlremmits grow by adding additional modxales. Refurbishment 
is Mcon^llshed by replacing modules, and development and qiuallflcatlon costs 
are zeduced. The radiator panels could en^iloy either flexible or xlgld fins, 
and heat pipe or punned fluid radiator panels could be used. The modules 
could be designed to fit into the Shuttle payload bay ^en zetracted, as shown 
In Figure 1$. For a 250 KW power module applicatlcxi, 8 to 10 modiiles could be 
Joined togeth^ as shown in Figure l6. 

3.2 Blanket Radiator System 

The blanket radiator system consists of submodtiles of heat pipes 
Joined by flexible radiator fin material as shown in Figures 17 and l6. The 
heat pipes may be Joined to rigid punned fluid manifolds by means of contact 
pressure to form large radiating areas as shown in Figure 19. This 
system has the advantages of modular designs described above, and may be 
stowed as a compact package. 

3.3 Elemental Heat Pipe Radiator 

The Idea behind the elemental heat pipe radiator, shown In Figure 
20 , is to build up spacecraft radiators from the simplest possible elements. 
The element consists of a single heat pipe connected to a radiator fin. For 
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snail heat rejection requirraente the heat plpe.i nay he tied directly to the 
heat source* vhereas for large heat reJectiMi requirenents* the heat pipes 
nay he connected to the heat source through a transport fluid loop as shown 
in Figure 21, The end of the heat pipe is specially formed to facilitate 
Joining the elenent to the heat source as is shown in Figure 22, for exaiiq>le. 
Series of eloaents nay he c-^nnected together as shown in Figures 23 and 2k 
to form large radiating surfaces. Typical design details for the elements are 
given in Figure 25, Advantages of this approach are listed in Table XVIII, 

3.1^ Extended Life Elexihle Radiator 

The extended life flexible rmdlator, shown In Figure 26, is a 
variation of the fl^clhle rmdlator previously developed under this contract, 
and could he huilt taslng existing technology. It achieves an extended operating 
life hy eaDploying metal tubing in place of the Teflon tubing of the soft tube 
flexible radiator. The deployment retraction system is also modified for 
long life. It retains many of the advantages of the soft tube flexible radian 
tor in that it may be developed and qualified Independent of its mission, can 
be stowed in a coa^act volume, is deployable and retractable, and has wide 
heat load capacity. It is compatible with most transport fluids, and has a 
broader operating toaperature range than the existing flexible radiator, 

3,5 Condensing Radiator 

In future applications where electrical power can be obtained for 
weight penalties on the order of 100 and mission durations are rela- 

tively long so that mlcrraneteorold protection for the radiator is a problem, 
or radiating surface area must be minimized, a condensing radiator system such 
as is lllxistrated in Figure 27 will have advantages over conventional radiators, 
Scxae of the operating features of the systma are listed in Figure 27, Vought 
has developed technology for fabricating this system under the Self Contcdned 
Heat Rejection Module program sponsored by NASA-JSC (MAS9-lUl»O0), 
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FIGURE 22 

ELEMENTAL HEAT PIPE RADIATOR MANIFOLD ATTACHMENT 
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TABLE XVI XI 

ELEMENTAL HEAT PIPE RADIATOR ADVANTAGES 


• DEVELOP AND QUALIFY INDEPENDENT OF MISSION 

• APPLICABLE TO LARGE OR SMALL SYSTEMS 

• STOW IN COMPACT VOLUME 

• LARGE SYSTEMS ASSEMBLED IN SPACE 

• REPAIR OR REPLACE WITHOUT OPENING PUMPED FLUID LOOP 

• LARGE NUMBER OF IDENTICAL ELEMENTS REDUCE COST 

• NO BRAZING OF HEAT PIPE REQUIRED 

• SIMPLE HEAT PIPE GEOMETRY 

• LOW THERMAL RESISTANCE AT JOINT 
ft. USES EXISTING TECHNOLOGY 
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• RADIATOR AREA REDUCED BY 50% 

• PUMP POWER IS 20% OF TOTAL HEAT LOAD 

• EXPANSION VALVE PROVIDES HEAT LOAD CONTROL 

• HEAT PIPE RADIATOR 

• CYCLES OTHER THAN VAPOR COMPRESSION NEED STUDY 


FIGURE 27 

CONDENSING RADIATOR SYSTEM 




4.0 LARGE LORO LIfE RADIATORS 

The optimization of large heat rejection systems involves the 
selection of components for reliability, weight, and performance; and the 
evaluation of edtextiative arrangements of these components with redundancy 
and/or oversizing to achieve a prescribed capacity and probability of mission 
success with minimum weigiit and cost. 

Because of the number of components required in large space heat 
rejection systems, and the failure rates of the components, redundant el^nents 
or subsystems are essential to hi^ system reliability. Figure 28 illustrates 
how redundancy ii^proves the probability of maintaining full operating capacity 
throughout a mission. In this example, each subsystem has a 90^ probability of 
s\tt*viving the mission. With one redundant subsystem (a toted of two subsystems) 
the probability that one or the other of the subsystems will fail is 18JI, and 
the probability that both systems will fail is 1?. Thus the probability that 
one of the two subsystems will remain functional is P = 1 - P(2) « 99^. For 
double redundemt systems the probability that one of the three loops will 
remain functional is P = 1-P(3) = 99.9)5. 

If, in the example with three subsystems, two of the systems must 
remain functional to maintain fvill system capacity, the probability of success 
is P = 1 - P(2) - P(3) “ 97*2)5. This illustrates how large systems can 

be constructed from modtales with oversizing to achieve high probabilities of 
success. 

There are several ways to provide redundancy. With full redundancy, 
each system element has a redundant element capable of performing the same 
function. With selective redvindancy, system elements are selectively grouped 
with a redundant identical group performing the same function. With simple 
redundancy, the basic system has a redundant system capable of performing the 
same function. In addition there are several ways to activate the redundant 
components. Table XIX gives advantages and disadvantages of alternative 
approaches. 

When constructing large heat rejection systems, there are cost and 
manufactxiring advantages in building up the system from a number of identical 
modules as described in Section 2.0. Hius the system would be based on simple 
redundancy although the modules themselves might employ selective redundancy 
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for coB®onent8 vlth poor reliabilities. Exsaole calculatione are given below 
to Illustrate how the system's probability of remaining functional varies with 
the nusiber of modules and the probability that each module will survive the 
mission. 

Figure 29 shows results obtained from analyses such as are illus- 
trated in Figure 28, for a system constructed from 10 independent modules. 

The figure shows the probability of the loss of subsystems as it affects the 
total radiating area of the system. The most probable percentage of area loss 
increases as the probability of survival of the subsystems decreases. If ihe 
system is oversised so that part of the radiating area can be lost without 
affecting the system's capacity to perform at the required level, the systesis 
survivability will exceed that of its subsystems. In this case the systems 
probability of stirvival is computed by summing the probabilities of all pos- 
sible combinations of sub^stem failures which will result in more area loss 
than is allowable. This is analogous to confuting the areas under the curves 
in Figure 29 past the point corresponding to the amount of system oversize. 

Figure 30 give? typical results confuted for a system with 20 subsystems. The 
figure shows that if a system survivability of 98 ^ is reqtiired (.02 probability 
of loss of more area than oversized), the system must be oversized by 9% if 
^he probability of stirvival of the subsystems is 0.99, by 13J( if the subsystem 
svirvivability is 0.95, and by 2b% if the subsystem survivability is 0.90. This 
means the system having 0.99 stirvivability modules must be designed to operate 
with l8 subsystems; the system having 0.95 modules with l6 subsystems; and 
the system having 0.90 modules with lU subsystems. Thus the size of the modules 
must increase as the probability of survival decreases. 

In some cases the edlowable radiating area oversizing may be limited. 
Thus Figure 31 shows how the probability of survival of a system with 20? 
oversizing depends on the number and survivability of its subsystems. This 
figure was generated from cross— plots of curves such as are given in Figure 30. 

The results show, for example, that a 98? probability of success can be obtained 
by a system of 20 modules, each having 95^ survivability; or by a system with 
10 modules, each having 98? survivability. Which system is better depends on 
properties such as weight, cost, and complexity. 

Figure 32 compares the weights computed for a large pumped fluid 
radiator system having variable numbers of subsystems and subsystem survivabilities. 
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FIGURE 30 EFFECT 








PROBABILITY OF LOSS OF MORE THAN 20% OF AREA 









The figure shove that the velgbt of the radiators decreases as the number of 
modules increases for a fixed probability of maintaining full i^stem capacity. 
Cost and other vehicle component weights mie^t increase for large numbers of 
modtiles, so that the system desiffi cannot be selected based on the radiator 
weight alone. However, the figure shows that signiflcent weight savings are 
possible if the syston is constructed from 10 or more modules. This number 
does not seem uxureasonable for a system of 200 KW capacity. 

General Conmarlsons of Pumped Fluid and Hybrid Heat Pi 


Fluid Radiators 

In this section the weights of radiator panels with and without 
heat pipes are cooqpared for various heat rejection loads and mission durations 
to determine whether the heat pipes reduce the radiator weight. Optimum designs 
of each type of radiator are compared. Several types of heat pipe radiator 
manifold/evaporator designs are considered to determine the effect of various 
degrees of added cong>lexlty on the toteLL radiator weight. Meteoroid bianper 
and meteoroid armor protection designs are clIso studied. In addition, alter- 
native ways of constzoicting largo systems from variable numbers of subsystems 
are considered for each type of design. Redundant flow passages attached to 
a common radiator fin are not considered. Redundant flow passages complicate 
the manifold/evaporator design of heat pipe based radiators so as to increase 
the temperature drop between the fluid loop and the radiating fin. Therefore, 
to avoid complicating the comparative evaluation of radiators with and without 
heat pipes, redundant flow loops are not considered. 

Redundant flow loops are considered in the more specific and de- 
tailed case study documented in the next section, where it is shown that the 
weight of punned fluid radiators can be reduced below the values reported here 
in some cases. The purpose of this section is to determine whether there are 
applications ndiere heat pipes offer advantages >dnen redundant flow loops on a 
common radiating surface cannot be considered. This a valid question for some 
potential future applications where redundant elements must be physically 
separated. Such requirements might be necessary to minimize the threat of 
damage to the heat rejection system by laser attacks, or gross wreckage by 
manipulators, etc. 



Figure 33 gives the radiator veldts* Including pua^lng power 
penalty for several desloas for a 20 KW systea with 99]( probability of not 
being incapacitated by micrometeorolds in a 9 year nlsslon. The heat pipe 
radiator entries are for the lightest weight of several heat pipe Manifold/ 
evaporator deelgns considered* This Is a proprletasy cracept generated under 
Vought funding. Other manifold/evaporator designs such as are common in 
the open literature are considerably heavier* 

The overall system probability of surviving 3 years In a micro- 
meteoroid environment is at least 99 % for each design represented in Figure 
33* For the designs with only one subsystem, the subsystem must have a 99 % 
probability of survival. With two subsystems, each having, a 90$ chance of 
surviving, the probability that one of the two subsystems will survive is 99$. 
The wall thicknesses of the 90$ subsystems are less than for the 99$ subsystems, 
but the wei^t of two 90$ subsystems is more than the veight of one 99$ sub- 
system, so that there is no advantage in such designs. The piunped fluid based 
designs with five subsystems, requiring three for full capacity, are weight 
competitive with the single subsystem designs, but the weight advantage is 
insufficient Justification for the added complexity. The designs with 10 or 
more subsystems have greater weight advantages, but probably would not be 
selected for a 20 KW system unless there were other reasons such as a need for 
redundant punq» loop components to improve the toted, heat rejection system 
reliability. 

Thus, for the 20 KW heat rejection requirement, the best system 
appears to be a single subsystem pumped fluid design with meteoroid bumpers. 

The beet heat pipe based radiator is lighter then the armor protected pumped 
fluid panel, but weighs more than the bumper protected design. Figure 3U shows 
radiating fin cross sections of* the three designs as determined frrai the com- 
puter routines described in Section 2. The weights shown in Figure 33 includes 
the wei^ts of the basic radiator panel components sho^m in Figure 3l» plus the 
wei^ts of the transport fluid, radiator coatings, manifolds, headers, pumping 
power penalty weight, and an additional structural weight. The additional 
structural weight accounts for honeycomb or other equivalent panel stiffening 
structure, euid is assumed to be Independent of the panel design* Figure 35 
gives the weight break-down for the three types of radiators. Figure ^6 gives 
radiator weights for a I60 KW system. Here the heat pipe based radiators weighs 
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FIGURE 34 20 KH RADIATOR DESIGN CONPARISONS. P(o) ■ 0.99 
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FIGURE 36 160 kW SYSTEM, 5 YEARS, P(o) « 


lesB than elthor the buiqper or antor protected designs unlese the ayetem Is 
built up from 5 or more subsystems. However, significant viifixt savings are 
possible if the system is constructed from 10 or more pumped fluid radiators. 
Similarly, the data for the 250 KW system in Figure 37 shove that large vei^t 
savings are possible if the system is constructed from subsystems of 10-20 KW 
capacity. 

Radiator veldts and siurface areas are plotted in Figures 38-It3 
for systems with 9^% probability of siurvlvlng five years. The results are 
similar to those for the 0.99 probability system discussed above excepb that 
there are no cases where heat pipe based desigas offer wel^t advantages. 

Thus, the results show that heat pipes axe best suited for large 
single subsystoB designs with high pzebabilities of surviving long periods 
of time on orbit (.99 for 5 years or more). For cases where the system can 
be built TQ> from modules of 10-20 KW capacity, piaiqped fluid designs axe 
lifter wei^t. 
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5.0 DESIOn OP A 50 Ktf HEAT REJECTION SISTEW 

In this section a complete heat rejection ayatem consisting of 
radiators and transport loop components is a.ialyzed to determine how the radia- 
tor design impacts the total system weight and reliability. The majority of the 
trade studies in this section were performed in conjunction with concurrent re- 
lated projects, and is summarized here to illustrate application of the techniques 
developed in the previous sections. The system is made up of two radiators (for 
symmetrical mass distribution) and associated hardware such as flexlines and fit- 
tings, fluid lines connecting the radiators to the heat load, pumps, motors, in- 
verters, check valves, accumulator, filter, fill and drain valves, temperature 
sensors, thermal control valves, and heat exchangers. Two fluid swivels are in- 
cluded for each radiator to permit continuous rotation of the radiators relative 
to the spacecraft. The system weight Includes a power peneJ.ty of 200 Ib/KW and 
a 0.65 Ib/ft penalty for structural support of the radiator. 

Figure shows typical radiator di men sions computed for this 
system. Since the radiator panel is relatively large, provisions are made for 
folding each panel in three places to facilitate stowage. Thus, with no 
redxindant loop, the heat pipe based panel shown in Figure would require 12 
littes (3 on each manifold). A pumped fluid radiator would require only 
two manifolds (or 6 flex lines) per panel. Extra manifolds are required in 
the heat pipe design because the lengths of the heat pipes cure constredned by 
watt-inch limitations. 

Table XX gives failure rate data for state-of-the-art spacecraft 
heat rejection system components. Table XXI illustrates how the failure rates 
data is used to determine the probability that the system will not fail in a 
5 year mission. For this example each Independent fluid loop ccn tains a re- 
dundant pump, and the radiacor is designed so that the probability of no 
meteoroid puncture Is 77^ for each loop. The probability that the system will 
remain operational ranges from 61tjt to 6d% depending on the number of flex lines 
required. Table XXII gives additional data showing how the system wei^t and 
reliability varies with the radiator construction, the probability of no 
puncture, and the number of independent flow loops. 

There are a few components in Table XXI whose relatively high 
failure rates add up so as to lower the reliability of the total system. 

Therefore, improvements in the system's probability of survival depend on 
whether these components can be Improved or backed up by redundant counterpeurts. 

The failure rate of radiators due to micrometeoroid Impact can be lowered from 
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^.6 • 10**^ to 0.22 * 10*^ Isy Increasing the wall tblcImesB of transport loop 
passages so that the probability of meteoroid penetration Increases from Tl% 
to 99JK. Similarly* the failure rate associated with meteoroid penetrations 
of connecting lines can be reduced significantly by Increasing the tubing wall 
thickness. Redundant thermal control valves and tes^erature sensors can also 
be provided without seriously Isqpactlng the cost or con^lexlty of the system. 
Complete redundant flow loops are required to avoid the limitations of flex 
lines or fluid swivels. 

Figure shows how the systems probability of failure can be 
reduced at the expense of added system coo^lexlty. Eadi of the flow loops for 
the systems represented In Figure contains redundant eluents for all com* 
ponents having high failure rates. Thus the system represented by the hexa- 
gonal symbols Is the some as that of Table XXII except that redundant thermal 
control valves and teiiQ>erature sensors have been added. The other systems 
have at least three Independent transport loops for added redundancy of cri- 
tical elements. 

Figure k3 shows that the heat pipe based radiators and pushed fluid 
radiators with meteoroid bunkers have practically the same wel^xt. However* 
the analyses assume that redundant; flow loops can be added to the heat pipe 
designs without Increasing the ten$>eratvire drop between the fluid loop and 
the heat pipes. It Is probable that uhen a design has been worked out to 
Incorporate redundant flow loops In heat pipe radiators that some Increase In 
radiating area and consequently weight* will be required to off-set the added 
thermal resistance between the transport fluid emd radiating fin. Thus the 
weights of the heat pipe based systems in Figure 4? are probably optimistic. 
Addltloncd costs for heat pipes and radiator panel fabrication wo\ild be 
another factor against the heat pipe design. Table XXIII summarizes some of 
the other factors that should be considered in trade studies for selecting a 
radiator design. 

In summary* the design analyses for a typical ^0 KH heat rejection 
system show that reliability problems inherent In pximp loop conqponents cure 
more serious than those associated with the radiator panels. Uhless redtmdant 
flow loops are employed* there is no Justification for designing the radiator 
with more than 90% probability of svirvlvlng In the mlcrometeorold environment. 
With redundant components the radiators should be teslgned for 95-^9% survi- 
vability. Flex lines and fluid swivels should be avoided whenever possible. 
There Is no apparent Justification for en5>loylng heat pipes In a ^0 KW radiator. 
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FIGURE 45 SYSTEM FAILURE PROBABILITIES AND HEIGHTS 

50 kH power module 
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X£F*XN/41.«'U3/04) ” 

F*4M: ♦XK3)/XEF 
IF4F1-F) 8U»70,7C 
7j F1«F 
60 TO S 
4)0 FlaF 

9'i 0XH*«0XM/2. 

IF4F2-FI 92i91»91 

91 F2»F 

60 TO 4 

92 F2 aF 

93 ODa-Ub/2. 

1F4F3-F) 95»9h»94 

94 F3>F 

60 TO 3 
9b F3aF 

96 OSs-US/2. 

1F4F4-F) 98,97*97 

97 FtaF 

60 TO 1 


5?'^ P«t » 


AS 




9d F«*sF 
99 aXM8-UXN/i. 

RReaf?E/l . 

XUsA/h/XL 

XSi»a<l/SF 

WRITttb»6c3l YL»<«X9.S»0.XL»XM,F 
MRltEltftt ) DIAN» JIAH.T,DELTAtOP»RRE«MT 
NRliEie>«b£C) idUt iC«T3*1*bH«« 8H«SH 
NRIlECutbh.) TtfS»S3tTS»0M 
MRlTE<t)»bbCI Wrt,WCV »HH,MF,WFL»WF 
•iRlTEIttObS) A,XEFfW,XLL»XSS 
100 CCNtINUE 
STOP 
ENO 


APPENDIX B 


LISTING OF COMPUTER PROGRAM FOR OPTIMIZATION 
OF PUMPED FUJID RADIATORS WITH BUMPER 
PROTECTION FOR MICIK)METEOROIDS 




PROGRAM OADOl lNPUT*OUrPUT« rAPB5»IMPUT,TAf*E6BOUTRUn 
OIMCNSION FUU*Q UJ 
54U format 3 > 

5j«» FO<MAr »F?. »,F6,1,4F6.3> 

90^ FORHUTun^) 

9u6 FO«‘1AT<iHl,ldH RADIAFOn F.UiO »»3Aiu»F) 

Slu FORHATd^icFe. J» 

601- FORMAri//f <»M s ,F7. lAh TIME «»EARS» s,Fr.3,7M P|m)s«F(>.3) 

61. FORk'AFU-H T IN ( F) = ,F 7. Et ilH F OUT(F) = ,F 7.i»li.h F AM'JtFI s«F7.&) 
6^, FOt<MAF^<♦^^ UELTA P <MAM IF OLP) 70FLFA P (CROSS TUBE! =«F6.3) 

630 FORMAT(7H P'JMP POWER PENALFt lES/M) a,F7.3l 
635 FORMAFUIH oumP EFF^SlENCV a,K6.3) 

64. FOPMAI(5Ji AREA FACTOR IN OPTIMIZATION PARAMETER (L8/FTZ) a,F8.4> 
650 F0RMAF(8H E 9ARa,F7.4»7H ZfcTAa»F7.4tl3M T COATING a, FA. 5) 

655 F0RMATI5N MJs,F»,j,. 5H LT/FT/HR RHOa,F8. 3,11H LB/FTi Ka,F8.4ti.5 
H D/^^/FT/• CPa,F7.4*fH 3 /lB/F) 

66v F0KMATI7(1«,F8,4)) 

66'i FORMAF(iXtr3.1«4 (iX « F8.4 ) « / I 

67w F0RMATC63H MAN OIA HEA3 OIA MALI. F FIN T DELTA P RE/lOOo 
• - 1 MF/A< 

68C FORHAFlSZM WT M MT CF HF H MT F MT FL MT PI 
69V F0RMAFI42H A EFF M NL NSI 

..-73u F0>HAF(4X»3M-/W*5X*5HN F IN«6X*lHSf 8X*1HD»8X*1HL»6X»5HM 00T.6X,1HP) 
REAn(5,505) FLUID 
REA3(5»514l XMU ♦ RH0,XK«CP 
- - READCS^Sj.I Q,FlME.Pfc 

tEAO(5t5:3) TiN»TUUT«TAMB 
RCAD(5*512) XK1*XK2«XK3 

REA0t5,57il XEfZ,TC ... 

READ(5«53jI XNP 
%EAO(5,510) NO 

Slua.l7l4/a0.**8l . 

PR* (CP*XMJ7if<xl^*. 33333 
RH0M=17 a» 

.... . XKNalJu. . 

PIa3* 141 59 
CMsRH0M*PI7l44. 

C(ial./Pl/X< 

CREa4l./Pl/XMU 

PLaALOG(P.| 

fI=TIN*46«. 

FOaFOUF*4bw. 

F AaTAMii^Aaj • 

01=ALOG( lfI-rA)*tTOflA)/(ri*TA)/cro-FA) ) 

DZi:AFAfl(TI7f A) - AT AN (TO/FA) 

DL=AwOGI (Fl*^4-TA**4</ (F0**4-TA**4I) 

04a(0*/4.-32/2.) /TA»*3 
CF»3.*SI..**t*( U rifT0)/2.) /XKM 

QaQ/i.. 

■ - DO 1C. Hgsl*N(] 

0=2. *0 

HRIT£(6,5.5» FLJIl) 

WRITE (b*695) X;iU» RHO« X<» 


6Z 




HtiTF (6f6u3 ) QtTIMEfPO 

riN»rouT,rAM.j 

KKi 

KRITH(6,63n 
NRITc UNP 

NI^irL owu ) x<» 

XF,Z,TC 

WRITi (u.Pju) 
write ( 6,6f J) 

MRITEI6,fiD;i) 

WKZrE(6,690) 

FSsl.J, 

XNL»L. 

IJi XNL=XRL*1, 

- • F 4 slii-. 

XNs.7 
OXNs.35 
00 99 1=1,2 
1 XN=XNfOXN 
CNsXN*XF*SI3 
Xa!l,-XN 

xs(,!)4t,z«^,Z9(:5*r^..ri.5*Y*r) *SQRr (tj 

5*2# 

03*1, 

FJaliia, 

DO 96 J*1,S 
■ 3 SaSfOS 

IF<S«6r,j«| 60 TO 1C3 
DS*0S/2, 

S*-OS 

U3 SF*S/l«:, 

OELTA=CF»<(SFTXM*t:,i 

- • . F2alCJ, . 

0 = 0 • , 6 
00a . .1 
00 93 K=l,i» 

4 OaQ^ao 

iF(0,6F,,,) 30 TO 1C4 
OUaOO/ 2 , 

Oa-03 

104 C0a*?i&/9H0/ ((1. ,*□) aaei 
XMaua.ao 

0XHax»4/2, 

Flai.3, 

OO 9w Lsl,5 

5 XHalMfOXi 

IF|XM,cr,3.) GO TO 1 j 5 
OXMaOXM/,;. 

XM=-OXM 

105 7 E = C.^E♦X^^/^ 

X NUs4 • 

XF*64*/R£ 


ORIIQIINAL PAGE IS 
OF POOR QUALITY 


0 


60 TO 

IFUr*6CJv«) 49«SJf4. 

3« yNUs.i.lb*(t?L***66ttT-i<j5«>*PR 

0Rel.«:?4 
60 TO 5, 

6U XNU«.4«.3*\F**.<J*Pr. 

XF*,*«<»/C<E**.2I 
- ■ 0Rs*.234 
51. RaCR/XNU 

03s5N*SP»<*D.. 

. CT»CU-T0l/<J«*^03i 

As5414.*0F0N/CT 
XL*XM»Ci»* tl I-TO) /CT/CN/SF 
... CAPtsXNL**„ 

«sA.'';apl 

TLHsCAPL/<« 

. .. DRa(DR*P**2.8/ <XXl*SF**i«d*XU)**.2w83 

OP*CO*XM*XM*X..»XF 

0IAMsOR*O 

IF(XNL-j«) 5t*Sl*5i 

-51 OIAHsOlAH - 

60 TO 53 

52 ORamCAPL-<U •♦1.8>*(CAPL»-4.6*XLI/XK1/M/(XL*»1*8I)**«2U8S 
0IAH».5146*0iAH*DR 

53 ARsn/SF»OIAM* { XNL».544I /CAPU2.*0IAH* ICAPL-XU FA 

... T®.l*«'»»<4-2*TlME»A»AR/PU**.2'a)/S5RTCRHOMl 

HP®.C3Ci;5423*<<2^XM*0P/RH0/SF/XL/XNP 
HHaCM»T*CT*(XNL*l.J ♦01AM*(XNL». 544H/CAPL 

riCTsCi4«T*(rf'))/SF 

MF=R40M* (TC»0FLTA)/12* 

RHaC<i*T*4 0IAH^Tl*2.*(l./W-XL/AI 

WFusniAM*3IAM» (i,/xL+.544FCAPU ♦0IAI-*01AH*2.* Cl*/M-XL/Ai ♦0»0/S 
RFL=*‘hO*PI»MFu/5T6. 

*(TsHM+wCT»rtH^WF4WFL->WP 

XEFaXN/U. ♦03/04) 

FsCHT^XKS) '<tF 
IFtFji-F) Tj,T.),7w 
-ifa Flap 
»0 TO 5 
Ou Flap 
Ota 0XMs>0XM/2* 

iFtFt-P) 9tt91t01 

91 F2»F 
60 TO 4 

92 FZaF 

9i DD=-DU/t.. 

IF4F3-F) 95,94,94 

94 F3*F 
60 TO 1 

95 F3=F 









OS«*OS/2. 

'1 60 TO 1 

’• 98 F4*F 

' 99 0XN«»3XN/2. 

RR6*RE/i'lJ3, 

’ Xi.L*A/H/X. 

XSSsH/SF 

MRirE«6«693l TUM,XN,S,0,XL,XH.F 
' HRXTE(6«86i<} 0IAM«0XAH*T«O£Lr4tOF«RFE«MT 

RRmC6,6o^) XI1,MCT,WH,WF,WFI.*WP 

WRITE I6«6S$) 4,XEF*««XLutXSS 
; IFIF5-F) iibtUhiUi 

" — a02.FS*F... - 

* 60 TO lOi 

1^16 F9«F 

■' XFIXNL*LT*4.5I 60 TO iJl .. 

10 a CONTINUE 
STOP 




« 
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IM 


_ me, r dT 
vour 




^ + la. 

X\ ^ TT 


4<r 




T’dr 1 

r*-T.‘) j 


i£T 


F = 


-T^ 


6 = 




i (T*~ T») 


H — F +■ ^ G vvhere d = Lt^ran^ multi pli 


EULEB5 EQUATION 

d- ( ^ = O 

dT I di i ax 














r 'w 


vj 


T^*dlT 




hJ(JM£R/CAL Solution 


jfi 

3,MPS0^S eoL£: |fcx)oU.4[fe.H4f0<,HK)] 


^tOO 

T dT 


- = T (<^-S5.S X lo'®) = 2.?% 


X 10 


«.5 


40 

/-/Oo 


s/ 


4o 


_ _L tf r S3S.704-1 
T*-T^ 4 . \_ 





" ).909C> X (o“^ 


C7 


s =rO,7 - 


(A 2351 > ~Uc ie \^I0)- 0,42857 


0.4Z6S1 

Vs/0)’ ■ ^ •^•^*** * 


AU’^viinAhk) 

Add) 


1 +d,25odKio' 

~Z.Z<»<, xio‘^1 

1"** 

|.44,lost 

JflO^fcX 10^ j 

1 +5.Z>o»a|^’^ 

o.nwsto 



<a 


